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Abstract 

I review the literature on observational aspects of waves in the solar chromosphere in the first part of this 
contribution. High-frequency waves are invoked to build elaborate cool-star chromosphere heating theories but 
have not been detected decisively so far, neither as magnetic modes in network elements nor as acoustic modes 
in below-the-canopy internetwork regions. Three-minute upward-propagating acoustic shocks are thoroughly 
established through numerical simulation as the cause of intermittent bright internetwork grains, but their pis- 
toning and their role in the low-chromosphere energy budget remain in debate. Three-minute wave interaction 
with magnetic canopies is a newer interest, presently progressing through numerical simulation. Three-minute 
umbral flashes and running penumbral waves seem a similar acoustic-shock phenomenon awaiting numerical 
simulation. The low-frequency network Doppler modulation remains enigmatic. 

In the second part, I address low-frequency ultraviolet brightness variations of the internetwork chromo- 
sphere in more detail. They contribute about half of the internetwork brightness modulation and presumably 
figure in cool-star basal flux. They appear to be a mixture of inverse-contrast granular overshoot at small scales 
and gravity-wave interference at mesogranular scales. I present TRACE evidence for the latter interpretation, 
and speculate that the low-frequency brightness minima map canopy heights. 

1 Introduction 

Waves in the solar chromosphere may be divided in many different ways: high versus low chromosphere, 
below the canopy versus above the canopy, network versus internetwork, acoustic versus magnetic versus 
gravity mode, linear versus shocks, quiet versus active region, standing versus propagating versus evanes- 
cent, global versus local, short-period versus three-minute versus five-minute versus long-period, longitudi- 
nal versus transverse versus torsional. And, in particular, whether they heat appreciably or negligibly. The 
latter issue is the most frequently quoted motivation for chromospheric wave studies and for their claims 
to fame, but it seems to me that identifying the nature of the observed modulations and the underlying 
processes and interactions should come first - and, personally, I regard the physics as more interesting than 
the eventually resulting overall energy budget, a mere summation at the very end with the answer already 
known as balancing photon losses. 



'Appeared in the proceedings of a NATO Advanced Research Workshop held at Budapest in 2002 edited by Erdel vTet all 1 20031) 
but not entered or even listed on ADS through laxness of the publisher and editors. Submitted to arXiv Astro-ph on December 6, 2010 
to remedy this neglect and to demonstrate my newest latex trick: in on-screen reading, clicking on the year in each citation opens the 
corresponding ADS abstract page in your browser. Explanation: |http: //www. astro .uu.nl/-rutten 
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2 Overview 

In this section I briefly review the literatur^U using the above splits as guideline. 



Chromosphere. Even this term i s no n-trivial. In the stand ard one-dimensional modeling of Vernazza, 
Avrett & Loeser dl973l lT976 19811) and Fontenla et a l. ( 1993), it denotes the plane-parallel layers between 
the temperature minimum and the onset of the coronal temperature rise. Within this definition, UV continua 
shortward of 1600 A come from the chromosphere (VALIII Fig. 36 in [Vernazza et alj 19811), and strong 
optical lines are also chromospheric - but not always. For example, the NalD lines reach r = 1 above 
the VALIII temperature minimum, but they are photospheric in their i ntensity response and do not map 
the chromospheric temperature rise in their VALIII emergent profiles dBruls et alj 19921). Even Call K 



filtergrams primarily sample the upper photosphere since the intensity-encoding thermal creation (rather 
than the velocity-encoding last scattering) of most of the observed photons takes place well below the 
VALIII minimum. In addition, the existence of a ubiquitous temperature starting at h — 500 km has been 



put into doubt, first by Ayres' CO modeling (e.g.jAvre s & Testerman|_l 98 II: lAyres 1981; Ayres et al. 1986; 
Ayres & Wiedemann! 1 9891: 1 Ayres & Brault 1990; Ayres & Rabinl 19961) and subsequently by the Carlsson- 
Stein acoustic shock modeling (see below). 

Hence, I prefer to define "chromosphere" - harking back to its eclipse origin = the purple color of Ha 
plus Hy6 off-limb emission - as the solar regime characteristically sampled by Balmer lines as a complex 
mass of fibrils, and taking the latter to represent (incomplete!) mappings of magnetic canopiefl 



Cano pies. The conc ept comes from IGiovanelli ( 1 9801) ! IJones & Giovanellil (19821) 



982), 
1993; 



Giovanelli & J ones ( 198 21) and is portraye d in older models and recent simulations (e.g jBiinte et al. 1993 ; 
Rosenthal et al.1 20021; iBogdan et a l. 2002) as smoothly upward spreading field hovering dome-like over 
essentially field-free internetwork. This is a simplification. If we define the chromosphere to start at canopy 
height (the height where the plasma beta drops through unity), its lower boundary will actually be a very 
warped surface, offset by dynamical flows and with large topological variations defined by the small-scale 
and large-scale field s trength and connectivity , as partially delineated by Ha as short local and long distant 
fibril connections (cf. lSchrijver & Tjtld 20021). 



High versus low chromosphere. This distinction now means well above and just above the canopy, 
respectively. Most of the references in this review pertain to the low chromosphere or even the up- 
per photosphere (say h = 400 - 800 km). Higher up, most of the recent wave literature em- 
ploys CDS and SUME R spectrometry to d i scuss to what height the internetwor k three-minute os- 
cillations penetrate (e . g., Steffens et a l. 1997; Carlsson e t al. 1997 ; ljudge et al. 1997; iDovle etaL 



Curdt & Heinze l 1998, GouttebrozeetalJ 1 999T[Doy le et aljl999i; |\Vikst0l et al. 2000; ljudge et al. 



Dovle et al. 


19981; 


Judge et al. 


20011; 



Mcintosh & Judge! 2001 ). Such penetration is likely to vary strongly with the actual field geometry. Oscil- 
lation analyses in tandem with moss and coronal loop diagnostics are yet scarce (but see Ineke de Moortel's 

contribution in this volume). 

Theoretical insight comes from beautiful Oslo simulations (iRosenthal et al. 2002: cf. 
Bogdan et al. 2002). I wonder whether such simulations for different stellar parameters may, at 



long last, explain the Wilso n-Bappu relation between stellar luminosity and Call H&K peak width 
(Wils on & Vainu Ba ppu 1957). It must describe fluxtube atmosphere properties, rather than non-magnetic 



'Note added December 5, 2010: All ADS abstracts for almost all solar physicists cited in this paper are collected at 
http : //www . astro . uu . nl/ ~rutten/solar_abstracts 

2 The definition implies that chromospheric studies must include Ho- even though its formation is singularly awkward through mix- 
ing thickness and thinness with NLTE opacity and source function sensitivities including Zanstra-like ionisation-plus-recombination 
photon conversions. An unpleasant but inescapable conclusion to one prefering clean lines such as Call H&K. Quantitative Hrr 
mapping requires considerable work on reliable Ho- interpretation, comparable to the ongoing efforts in Stokes profile inversion. 
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atmospheric stratifications (e.g., lAvresI 19791; |Kneerj 19831), since the emission peaks come from the 
magnetic component. 



Network versus internetwork. This division is also overly simplistic. The verdict on internet- 
work fields isn't yet in, but they exist undoubtedly at some level of field strength and scale 
of spatial organization. Call K and TRACE UV image sequences show extended zones of en- 
hanced brightness (with respect to the darkest "cell centers") as aureoles around network (cf. 
Fig. [5}. Acoustic maps show similar thr ee-minute power aureoles when sampling power below the 
canopy and power shadows above it (e.g., iBraun et al. 1992; [Brown et al. 1992; iToner & Labo nte 1993; 
Hindman & Brownl 19981; iLindse v & Braun 1999; Braun & Lindsevf 19991; iThomas & Stanchfi eld 2000; 



Mcintosh e t al. 2001 [judge et alj 200li; lMcIntosh & Ju dge 2001 ; Kjij£gretalj200li) 



In addition, Ca II K and TRACE UV image sequences also show evidence of "magnetic fl ashers", pre- 
sumably isolated fluxtubes on their way to or from the ne twork concentrations (Brandt et al., 1992 1994; 
iNindos & Zirir] 1998l; lLiteset al. 1999; KjijgeretalJ 20011). Their isolation may be helpful in trying to iden- 
tify chromospheric tube modes excited by convective buffeting without amplitude loss from phase mixing 
over multiple fluxtubes. 



Acoustic versus magnetic versus gravity modes. Or mixtures, of course. Magnetic modes 
are often invoked to c o nvey energy from photospheric buffeting up along network fluxtubes (e.g., 
Hasan & Kalkofen 1999; Hasan et al. 2000) but have not been convincingly diagnosed yet, except 



negatively as suppression of acoustic wave power by conversion at canopies dMcIntosh et al. 2001 ; 
Mcintosh & Judgef 200 lj; [Rosenthal et ail 20021). 

Interna l gravity waves should be "c opiously excit ed " in granu la r overshoot acco rding to the- 
ory (e.g.. IWhitakel[T963i Lighfhill in iThomasI 19671; 
Toomre 19811 19821), but they are difficult to detect. 
edly. 



SteirJ 19671; ISchmiederi 19771; Mihalas & 
propagating slant- 



being sma ll -scale and 

The obse rvation al evidence is mostly indirect (jFrazierj 1968; Schmieder 1976; Cram 1978 



Brown & Harrison 1980; Durrant & Nesis 1981 ; IStaiger e t al. 1984; IStaigeri 19871; iDeubner & Fl eck 1989; 
Bonet et al. 19911; Komm et al. 19911), with the cl earest demonstration coming from wavenumber- and 



frequency-resolved {kh,f) phase-difference spectra (IKneer & von Uexkull 1993; Stra us & Bonaccinil 19971; 
Kriig er et alJ 2001 ). Whether they affect upper-photosphere or low-chromosphere energy balances is not 



known. See also Section[3]below. 



Linear versus shock behavio r. At lea st in the case of the three-minute oscillation, protests 
against large non-linearity dDeubnerl 19911) were silenced by the successful repro ductio n of so - called 
Call H?y grain behav ior in the celebrated simulations of Carlsson and Stein (1994| 119961 1997; 
Stein & Carlsson 1997) of the spectral time sequences of iLites et al. 1(19931) The Carlsson-Stein re- 



production of complex spectral Call H core evolution patterns identified the three-minute oscil latio n 
beyo nd doubt with upward pro pagating acoustic shock t rains, as proposed ea rlier by Athay (1970) , 

\Crcm ( 1 Qll\ T in Mr QViimanir-Vil ( 1 Q7Zlh llVToin »t al I M Q871 T olKa^lior at al I (1 Button St I Ti tonKi-rvcO- 



ICraml dT972l)l lliu & Skumanichl ( 19741) Mein e t al. ( 19871) iLeibacher et al.1 ( 19821) Rutten & Uitenbroek 
( 1991b 1991a) and Rammacher & Ulmschneider ( 1992) However, the nature of the photospheric piston- 



ing remains in debate dN indos & Zirin 19981; ILites et al. 1999 ; IWorden et all 19991; ISivaraman et al.1 2000 ; 
Skar tlien et all 2000 ; Hoekzema et al. 2002 as does the role or absence of a role in c hromospheric heating 
dCarlsson & Steinl 1995L iTheureret all 1997al; iKalkofen et ail 19991; iKalkofenl 20011). See also Section |3] 
below. 



Active region oscillations. A direct active-region counterpart to the quiet-sun three-min ute oscillation 
consists of the three-minute oscillations produc ing umbral flashes (Beckers & Tallantl 19691). They show 
even more outspoken nonlinear character (e.g., Thomas et al. 19841; iLitesI 19861). and are most likely sim- 
ilar nearly-acoustic shock trains running up along the radial fields above umbrae as suggested by Lites 
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(1992 19941) So me are observed to p enetrate up into coronal loops in U V and EUV data (e.g., Brynild- 
sen et al., 1999 200 21; Maltby et aL, 199 9! 20011; JSchriiver et al.1 19991) and in microwave observations 

dGelfreikh et alj 1999l; ISMbasakil 2001l; lNindos et alj 20 021). _ _ 

Similarly, run ning penumbral waves (Zirin & Stein 1972; lAlissandrakis et al ] 1992; 
Brisken & Zirin 1997; Christopoulou et al., 2 0001, 20001, 20011, 2002) may harbor complex wave re- 



flection and mode conversion as postulated by Bogdan et al. 2002 from comparison to thinner fiuxtubes, 
but recent high-resolution data from three solar telescopes on La Palma suggest strongly that running 
penumbral waves and umbral flashes are both upward-propagating acoustic shock trains, seen differently 
through difference in field alignment with the line of sight and with the ap parent horizontal spreading 
pr imarily a mapping of the field geometry (IRouppe van der Voort et al.1 20031). Oslo simulations as those 
of lCarlsson & Steinl ( 1997 1) are needed for dec i sive d iagnosis, including "calibration" of the polarimetric 
umbral-flash inversions of iSocas-Navarro et al.l (2001) through less indirect forward modeling. 

T he existence and nature of oscill at ions in the su nspot magnetic fi eld itsel f remain contested 
(e.g., iLites etalJ 19981: iBalfhasari 19991; [staude 20021: iKupke et al.l 2000). with ISettele et al.l ( 20021) 



warning against instrumental crosstalk and initial results fro m infrared spectropolarimetry indicat - 
ing that opacity variati ons are an alternative explanation (cf. Collados 2002; iKhomenko et al. 20031). 



Georgakilas et al. (2002) found indications from TRACE data that running penumbral waves even spread 



into sunspot moats. 

Chromospheric oscillations above plage have not received much attention in the literature, but Rita 
Ryutova's presentation of beautiful TRACE 171 A time slices from Richard Shine show interesting braiding 
with shorter braid period at larger magnetic filling factor (see her contribution in these proceedings). 



Standing versus propagatio n. This issue is also debated in the three-minute o scillation literature (Deub- 
neret al., 1992 19961: ISteffens et al.1 1995l: lHofmann et al.l 1996l; lDeubner 1998 ). with simulations showing 
that at least some of the observed standing-wave behavior comes from reflection s off the grain-forming 
shocks themselves and off canopies (ICarlsson & Steinl 1999l; lRosenthal et al. 2002). 

Another standing versus propagation issue concerns the nature of the so-called pseu do-ridges above 
the cutoff frequency in (ky,,f) diagrams, also a global-versus-local issue. Kumar (19941) described them 
as interference between directly emitted and once-bounced outgoing waves (that propagate up rather than 
being evanescent), using "interference" as a misleading mathematical term describing low-/ Fourier de- 
composition. Physically there is no actual wave interference, but simply power addition in (£/,,/) dia- 
grams at those (fe/,,/) locations that sample one-bounce horizontal spatial wavelengths at the correspond- 
ing temporal fr equency, just anot her expression of the single-bounce three-min ute power ridge in time- 
distance plots dDuvall et all 19931) and exhibiting the Duvall dispersion law dDuvalll 1982; Eq. 2.14 of 
Deubner & Goughl 1984). Below the cutoff frequency, the acoustic ridges describ e evanescent p-mod es, 



with small phase delays compared to the photosphere governed by non-adiabaticity dKrijger et al .1 2001 



Short-period versus three-minute versus five-minute versus low-frequency modulation. These terms 
are all often misnomers. By short-period waves one simply implies those that might after all heat the chro- 
mosphere, with Peter Ulmschneider with coworkers as tenacious champion. After giving up on five-minute 
heating when the /?-modes were identified and on three-minute heating when the Carlsson-Stein simu- 
lation refuted a low-lying chromospheric temperatu re rise, his quest turned to hi gher-frequency compo- 
nents not present in the Carl sson-Stein piston (e.g., Sutmann & Ulmschneidei 19 951; Theurer et a l., 1997a 
1997bl; iKalkofen et al. 1999) with obvious stellar overtones dUlmschneider et al.l 19961 ICuntz et al.[j 999). 
and m ore recently to adding longitudi nal, transver s e and t orsional tube waves (Ulmschneider et al., 2001b 
200 lal, Musielak and Ulmschneider, 2002, 2002a, 2002b) as well as detailed theoretical prediction of the 
mix of acoustic and tube waves that should explain observed cool-s tar chromospheric pho ton losses with 
the magnetic filling factor as activity scaling parameter (Fawzy et al., 2002a, 2002c, 2002b). 

Heating by 10-100 mHz waves, outside and inside fiuxtubes, is probable; the question is to 
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what extent. Observationally, they are hard to see since they require fast cadence and high angu- 
lar resoluti on and suffer response function loss b y vertical wavelengths fitting in contribution func- 
tions (e.g.. [Sc hmieder & M eirl 19871 lDurranj | [T980i) and thro ugh large sensitivity to seeing noise (e.g., 



Endler & Deub ner 1983; IDeubner et alj 19841: lutes et al.l 19941). Indeed, the tri al s to detect high-frequency 
waves remain inconclusive (e.g., "IPe ubner 19761 iLites & ChipmaijfffTl iMein & SchrruedeilT98TI 



Deubner & Fleckl 1990; Krijger_etaL|2001). 



The chromospheric three-minute and five-minute oscillations are discussed above. The names imply 
broad-band frequency domains, not specific periodicities but just a shift in dominance from five-minute to 
three-minute periodicity when rising from the photosphere up to the chromosphere (in the internetwork at 
least). It is therefore also wrong to call them chrom ospheric and photos pheric, respe ctively. The original 
Leighton-Simon-Noyes era diagram in Noyes ( 1967 ), reprinted in Rutten ( 1995 , 2001 ) clearly demonstrates 
the gradual change. 

Finally, low-frequency oscillations is a misnomer when the observed modulation is not oscillatory. This 
may hold for the observed low-frequency Dopplershift power of c homospheric network which may re- 



sult from convective fluxtube buff eting (cf. Kneer and von Uexkiill 1985 1986: 1 von Uexkiill e t al. 1989 



Lites et alj 19931; lHasan et al . 2000), and also for low-frequency low-chromosphere internetwork oscilla- 



tions if granular overshoot rather than internal gravity waves dominates the observed low-frequency power. 



3 Low-frequency internetwork background modulation 

I devote the remainder of this contribution to the topic on which I concentrated in my oral presentation at 
the meeting: the nature of the low-frequency background in low-chromosphere ("upper photosphere" in 
the under-the-canopy definition) image sequences. I have long been puzzled b y the mesh-like background 



pattern underlying ¥±2v grains in Ca II K image sequences (cf. lLites et alj 1999 ). The upshot is that I believe 
the answer to be a mixture of granular overshoot plus gravity-wave interference at slightly larger scales. I 
advoc ate this i nterpretation with selected results from TRACE data analyzed by J.M. Krijger in his PhD 
thesis (Krijger 2002 ). The diagrams in this section als o illustrate various p oints made abovqj. 

The observation and reduction details are given in lKrijger et alj (20011) The displays below come from 
the very quiet-sun data taken on October 14, 1998, when TRACE registered image sequences in its 1700, 
1600 and 1550 A ultraviolet passbands and also in white light. The combination permits comparison of pho- 
tospheric brightness patterns to the co-spatial and co-temporal ultraviolet ones. Such comparisons used to be 
made with groundbased telescopes combining Call K and continuum image registration (e.g., Hoekzema 
et al., 1998 1998 2002) but TRACE furnishes better quality thanks to the absence of seeing in spacfl 



Space-time representations. Figure Q] defines the topic. It compares the white-light low-photosphere 
scene with the co-spatial and co-temporal 1700 A high-photosphere scene in the upper panels. The lower 
panels compare evolutionary characteristics between the two regimes. The displays sample only very small 
parts of the full white-light and 1700 A data cube^fl 

The difference between dark internetwork and bright network is obvious in the ultraviolet image (sec- 
ond column). The network grains stand out even though their brightness is cut in half by the display 
scaling. TRACE'S resolution (1 arcsec) is insufficient to resolve the corresponding white-light "network 
bright points" (magnetic elements) residing within the underlying intergranular lanes. 

3 Note added on December 5, 2010: part of this work was published in Rutten & Kriiger (2003) but not all diagrams and ideas 
given here survived the referee. 

4 But the tide turns, now that speckle reconstruction, phase-diverse restoration and adaptive optics also do away with (most of) 
the seeing, permitting higher angular resolution than TRACE'S 1 arcsec. Much sharper images now result e.g., from our Dutch Open 
Telescope which is presently being equipped with multi- wavelength tomography capability I http : //dot . astro . uu . nil. 

s These cutouts are small in order to obtain sufficient magnification. They represent a very limited rendering of the 
full data cubes. A better view of the dynamical ultraviolet scene is gained by inspecting the TRACE movies available at 
http: //www. astro.uu.nl/-rutten/tracel 
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Figure 1: TRACE image and time-slice samples. First column: white-light intensity. Second column: 
co-temporal and co-spatial 1700 A intensity. The images in the upper row are small cut-outs of the full 
256x256 arcsec 2 field shown inFig.\5\ The time of observation is halfway the time slices in the lower row, as 
indicated by white markers. The slices show the temporal intensity evolution during 30 min of observation, 
for a horizontal cut through an internetwork area indicated by the white markers in the upper panels. It 
passes through weak network at the right. In the first two columns the greyscale is linear for all four panels, 
but it has been clipped at half the actual maximum for the 1700 A image to enhance internetwork. Third 
and fourth column: the same, but low-pass filtered (subsonic horizontal propagation only). Third column: 
white-light brightness on a sign-reversed greyscale. Fourth column: logarithm of the 1700 A brughtness 
temperature. 

The white-light time slice in the lower-left panel displays primarily granular evolution, with some larger- 
scale five-minute /?-mode modulation. The ultraviolet time slice shows the dynamical behavior characteris- 
tic of internetwork in the form of ubiquitous short-lived three-minute oscillation sequences. The brightest 
phases are called internetwork grains. The weak network grain near X = -100 arcsec stands out by its 
relative longevity. Note that the cut selection favours internetwork; stronger network grains produce much 
brighter vertical streaks. 

The two righthand columns show the same data after Fourier filtering and with modified greyscaling. 
In the third column (white light) low-pass "subsonic' filtering, i.e., applying a 3D Fourier "cone" filter to 
the transformed data cube which passes all signals with apparent horizontal speed below the sound speed, 
has removed the photospheric five-minute oscillation so that only the granular evolution patterns remain. 
In addition, the greyscale is sign-reversed to simulate "reversed granulation" with reversed contrast. The 
resulting mesh pattern in the reversed image illustrates the topological difference between granules and 
intergranular lanes. The corresponding time slice shows intergranular lanes as rather long-lived bright 
streaks. 

In the fourth column (1700 A) the low -pass filtering has removed the chromospheric three-minute os- 
cillation and therefore emphasizes the slower background evolution. The background streaks are relatively 
short and show larger horizontal displacements (tilts) than the intergranular streaks. At this resolution, there 
is no obvious correspondence between the reversed low-frequency granular evolution pattern in the third 
slice and the low-frequency internetwork background modulation pattern in the fourth slice. 
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Figure 2: Strous-fonnat scatter diagrams. The crowded central parts are plotted as sample density contours 
to avoid plot saturation. First row: time-delayed white-light brightness temperature in K WLit + At) against 
WL(t). Second row: time-delayed 1700 A brightness temperature UV(t + At) against UV(t). Third row: 
UV(t + At) against WL(t). The time delays At are specified in each panel and increase from left to right. 
The solid curves in the first panels show the occurrence distributions of the quantities plotted along x and 
y on inverted normalized scales. The dashed curves show the first moments of the sample density along 
horizontal and vertical cuts through the contours. 



Time-delay scatter representat ions. Figu res [2|-[3] are pixel-by-pixel dual-image correlation plots in an 
informative format initiated bv lStrousI ( 19941) Each plot is measured from a 256 x 256 px 2 TRACE subfield 
after removal of Fourier taper edges, image conversion into brightness temperature, and 1 .5 x 1 .5 arcsec 2 
boxcar smoothing to suppress noise. 

For each pixel in the subfield, i.e., each solar location, the brightness temperature at one moment in the 
one type of image (say white light = WL) is taken as x quantity plotted horizontally, its value at another 
(later) moment in the other type of image (say 1700 A = UV) as y quantity plotted vertically. The pair 
defines one point in a scatter plot. Such pairwise comparisons are made for all pixels (solar locations) in 
the whole subfield or in some selective part of it, and repeated for 50 consecutive image pairs to gain high 
significance (millions of spatio-temporal samples). Plot saturation (total blackness) is avoided by plotting 
sample density contours instead of the individual pixel-by-pixel samples. 

The WL-WL comparison in the first row of Fig. [2] illustrates the format. Brightness distribution curves 
are added in the first panel. They are virtually the same for the other panels and also along x and y in 
this auto-correlation sequence. For 100% correlation all samples and both first-moment (center of gravity) 
curves lie along the forward diagonal. For 100% anticorrelation they lie along the backward diagonal. 
At the absence of any correlation the first-moment curves become perpendicular, parallel to the axes, and 




Figure 3: Time-delay scatter diagrams as in Fig. \2\butfor internetwork (IN) only. First row: time-delayed 
low-pass filtered 1700 A brightness temperature UV(t + Af) against low-pass filtered white-light brightness 
temperature WL(t). Second row: unfiltered UV against low-pass UV. Third row: low-pass UV against 
low-pass UV. Fourth row: high-pass UV against low-pass UV, with the mean UV internetwork value (T\, = 
4418 K) added to the high-pass UV modulation for y-axis compatibility. 



the contours become circular if the distribution function is symmetric. From left to right the panels have 
increasing time delay between the sampling of each type of image per solar location. The initial panel for 
Af = 22 s shows very high pattern correlation because the granulation has not changed much during this 
brief interval. The final panel for Af = 30 min shows absence of pattern correlation in the form of non- 
aligned first-moment curves and a nearly round bull's-eye contour pattern. The sequence illustrates that 
granulation largely loses its pattern identity over ten minutes, completely in half an hour. 

The UV-UV autocorrelation sequence in the second row illustrates the two-component dichotomy be- 
tween chromospheric network and internetwork. The bright upward distribution tail made up by network 
grains persists over long delays, whereas the darker internetwork (lower-left contour mountain) shows faster 
pattern change. 

The UV-WL cross-correlation sequence in the bottom row of Fig. [2] shows some persistent bright- 
bright correlation due to network, slight anticorrelation for low WL at Af = 3 min, and subsequent lack of 
persistent correlation at low UV brightness. 
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In Fig.[3]the ultraviolet signal is decomposed in constituents by using selective data subsets for similar 
scatter diagrams. All four rows are limited to internetwork (IN) areas only. Fourier cone filtering as in Fig.Q] 
is applied to show low-pass signals (apparent horizontal propagation below 6 km s -1 ) or high-pass signals 
(above 7 km s _1 ) only, selecting acoustic and non-acoustic modulation. The corresponding reduction in 
sampling statistics is offset by extending the pixel-by-pixel sampling to 120 consecutive image pairs per 
comparison (in all scatter plots the outer contour lies at 100 samples per bin with 25 bins per axis). The 
outer scatter clouds are now excluded to reduce the plot file size. 

The first row of Fig.[3]plots low -pass UV against low -pass WL. There is a slight but persistent anticorre- 
lation over time delay Af = 2 - 6 min. I attribute it to reversed granulation and expect it to become stronger 
in higher-resolution data that resolve granulation better than TRACE does. The slightness of the correlation 
shows that at the somewhat larger meso-scales imaged properly by TRACE, something else is present that 
does not obey point-to-point correlation or counter-correlation with the underlying photosphere at any time 
delay. 

The remaining three rows decompose the UV signal in high-frequency and low-frequency components, 
as contributed by the latter (low-pass UV as x-axis quantity). The second row plots total UV against 
low-pass UV in internetwork. The high initial correlation, with a roughly 2:1 slope, shows that peaks in 
internetwork UV brightness, i.e., internetwork grains, occur preferentially when a lso the slow-changing 
internetwork background is bright. This strong correlation is conform the finding of ICram & Dam e ( 1983) 
that bright Ca II H2V grains are invariably part of a larger-scale modulation pattern. The same happens 
in the ultraviolet continua. Thus, the slowly evolving background pattern and the three-minute oscillation 
combine constructively to produce bright internetwork grains, each contributing about half of the excess 
brightness tempera ture. This makes the b ackground a much more important grain co-localize^ than e.g., 
acoustic events (cf. iHoekzema et al.l 20021). 

The next row shows a low-pass internetwork-only UV autocorrelation sequence. It illustrates the five- 
minute lifetime of the low-frequency pattern at a given spatial position. Horizontal propagation - seen as 
tilts in the low-pass UV time slice in the final panel of Fig. Q] when having an x component - causes loss of 
scatter correlation when the travel exceeds the 1.5 arcsec boxcar smoothing. 

The final row correlates high-pass UV with low-pass UV in internetwork. The curved shape of the 
vertical moment curves indicates that large three-minute excursions, both to high and to low UV brightness 
temperature with respect to the mean after removal of all slow variations, tend partially toward correlation 
with large UV background amplitude. This implies that the three-minute oscillation partially has modula- 
tory character, gaining larger oscillation amplitude at larger background amplitude. The effect is not large. 
In addition, the contours have slightly asymmetric shape indicating non-linearity in wave behavior. When 
plotted as intensity, the scatter clouds stretch much further upwards due to the non-linear Planck function 
response in the ultraviolet. The conversion to brightness temperature circularizes the contours considerably. 

The conclusions from Figs. |2]-[3]are, first, that the slowly-evolving internetwork ultraviolet background 
contributes about half of internetwork grain brightness excesses, in addition to the acoustic modulation 
modeled by Carlsson & Stein, and, second, that at TRACE'S meso-granular resolution this background 
seems to be something else than reversed granulation. 



Fourier representations. Figure |4] compares the same TRACE data in terms of (£/,,/) Fourier phase 
differences, at left between 1700 and 1600 A brightness, at right between white-light and 1700 A brightness. 
The acoustic oscillations (ridges and pseudo-ridges) in a similar diagram from the May 12, 1998 TRACE 
data are discussed in Sect. 4.4 of Krii g er et al.l ( 200 II) ; the emphasis here lies on the low-frequency domain. 
Note that the October 14, 1998 TRACE data contain only very weak network (cf. Fig.|5]l so that these (£/,, f) 
diagrams are dominated by internetwork behavior. 

6 The low-frequency background was not part of the Carlsson-Stein simulation. The corresponding Dopplershifts have less power 
than the brightnes s modulation and were not part of the Carlsson-Stein iron-blend piston because they were removed in spectrograph 
drift correction by Lites e t all 1 19931) 
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Figure 4: Two-dimensional Fourier spectra from TRACE. Left hand graph: A0(17OO-16OO) intensity phase 
difference. Axes: horizontal wavenumber kh and temporal frequency f. The corresponding wavelengths and 
periodicities are specified along the top and righthand side. Greyscale: phase difference coding specified 
in the bar above the graph. The white curves along the sides are the temporal and spatial means, on linear 
scales in arbitrary units. The contours specify 1700- 1600 A coherence at values C — 0.4, 0.6, 0.8 and 
0.95, with the ticks directed to lower values. Righthand graph: similar A0(WL — 1700) phase difference 
spectrum between white light intensity and 1700 A intensity. The contours specify coherence levels C —0.1, 
0.2 and 0.5. Thijs Krijger production. 



The lefthand graph shows a prominent wedge-shaped signature of negative propagation over a large kh 
range. It has high coherence between 1700 and 1600 A modulation. The wedge location, its Lamb-line 
delimitation, and its negative values all suggest atmospheric gravity waves as cause. The large coherence 
implies that these waves dominate the low-frequency internetwork background in the ultraviolet at meso- 
scale wavelengths. 

The righthand (kh, f) diagram is very noisy but again shows a low-/ low-A; wedge of negative phase for 
kh < 2 arcsec -1 , i.e., at mesogranular scales. It reaches much larger A(p values due to the much larger span 
in formation height and it has much smaller coherence, but it is qualitatively similar. At larger kh the low-/ 
phase differences flips from large negative to large positive values at about kh = 2.5 arcsec -1 , attributed to 
wraparound in the A(f> = [-180, +180] deg evaluation. This makes the positive blob around kh — 3 arcsec" 1 
(outlined by a C = 0.1 coherence contour) a continuation of the negative dark wedge, presumably marking 
reverse granulation. The slight counter-correlation signature in the top row of Fig. [3] is therefore indeed 
caused at granular scales. This high-A: contribution is likely to become better defined at higher angular 
resolution than TRACE. The mesoscale negative-phase blob in the righthand (kh, f) diagram describes the 
non-correlated other agent, which I attribute to internal gravity waves. 



Discussion. In summary, it seems that the slowly-evolving background mesh pattern in the internetwork 
parts of chromospheric image sequences is made up of reversed granulation and internal gravity waves at 
slightly larger scales. 

Numerical simulation may address this interpretation beyond speculation. Ob viously, this should be fea- 
sible with realistic numerical simulations of the s olar granulation such as those of lStein & Nordlundl (1998) 
when equipped with sufficient atmosphere as in lSkartlien et al.l (20001) They should reproduce both granu- 
lar overshoot and gravity-wave emission. The two phenomena are akin but not identical, since overshoot is 
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a local phenom enon whereas gravity wave s spread away from their source, interfere, and break at relatively 



large height (cf. Mih alas & Toomrd 19811) 



Basal flux. The demonstration in Fig.[3]that roughly half of the brightness excess in internetwork grains 
is contributed by the low-frequency variation applies likewise to Ca II K.2V grains and the accompanying 
Ca II K line-core and inner- wing brightness modulation. Since the solar internetwork scene observed in 
Ca II K is dominated by the superposition of these acoustic and low-frequency variations, the corollary is 
that part of the basal flux observed from non-active cool stars Schrijver 1995, and references thereml is 



attributable to similar combination of acoustic and low-frequency phenomena including gravity waves - 
not just acoustics alone. 

Internetwork fields. As noted in Section|2]the debate on the pistoning of internetwork grains continues. 
It includes pro-and-contra claims on internetwork fields as grain localiser. Advocates pro may argue that 
the 50% grain co-localization by the low-frequency chromospheric background pattern seen in Fig. [3] re- 
sults from weak fields that are swept into in tergranular lanes wit hout reaching the magnetic flux of network 



elements, as sugge sted observationally bvlLin & Rimmele ( 1999) and theoretically by the magnetoconvec 



tion simulations of lEmonet & Cattaneo (2001), on the assumption that such small concentrations share the 
network habits of being bright in the chromosphere and displaying low-frequency variations there. 

I rather doubt this alternative explanation because of the ubiquity and regularity of the background 
emission, the low counter-correlation to intergranular lanes for the spatial "meso"-scales producing the 
negative-phase blob in Fig. [4] and the issue how intrinsically weak fields would influence the thermody- 
namics of the low chromosphere. The inverse-amplitude mapping in Fig.|5]discussed below suggests rather 
the reverse to me, namely that larger low-frequency modulation amplitude implies weaker magnetism. 

Nevertheless, internetwork magnetism cannot be rejected as a chromospheric brightness producer. 
I wonder whether it may play a role in the generation of small-scale reversed granulation, in the 
form of kilogauss fluxtubes that are too thin to be magnetographically observable so far. One rea- 
son for such wondering is the wide extent of relat ively bright areas surrounding chromospheric network 



the "intermediate" pixel class of lKriiger et al.l 20011 suggesting searches for relations between the occur- 



rence of reversed granulation and magnetic flux density in high-resolution large-sensitivity magnetograms. 

Inverse canopy mapping. Finally, Fig. [5] serves to illustrate a speculation bringing me back to the first 
item in my review above: the definition of the canopy as lower chromosphere boundary. It shows a one-hour 
temporal average of the 1 700 A brightness over the full TRACE field, with inverted ( 1/value) grey scale. The 
inversion de-emphasizes network grains by making them dark. The hour-long averaging reduces the three- 
minute component; remaining enhancements from internetwork grain nonlinearity (non-sinusoidal spikes) 
also become dark. Thus, the brightest features in this display mark locations where the low-frequency 
background modulation reaches its deepest minima without contamination by other phenomena. 

Figure [5] shows a grainy pattern. The brightest grains mark preponderance of extreme low-frequency 
minima. In reversed-granulation terms, these should correspond to long-lived or repetitive bright granules 
surviving a full hour of temporal averaging. Their distribution over the many network cells covered by 
Fig- clearly shows preference for the internetwork centers, with striking avoidance of the cell boundaries. 

My speculation is that such preferential occurrence of low-frequency extrema marks locations where 
the magnetic canopy reaches larger height than elsewhere. In the gravity-wave interpretation, amplitude 
reduction may occur where the field becomes strong enough to convert gravity waves into other modes or to 
reflect them. The areas containing the brightest inverse grains in Fig.|5]would then have the highest canopy, 
giving the waves more room to increase their amplitude - I suspect up to wave-breaking height. 

A reverse speculation by our meeting's co-director R. Erdelyi is that magnetic resonant absorption 
enhances spec ific wave modes at the canopy hei ght and that large apparent amplitude may mark such 



resonance (cf. lCadez et alJ 1997l; lPinter et al.l 2001 ) 
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Figure 5: Temporally-averaged 1700 A brightness on an inverted grey scale. The white lines mark division 
into smaller quadrants for processing. The yet smaller subfield of Fig. \I]is near the lower-left corner. The 
horizontal striping near the right edge results from solar rotation correction. 



4 Conclusion 



Both the review in Section [2] and the ultraviolet background analysis in Section [3] are observation- 
ally oriented. However, the major advance in chromospheric wave research in the past years is, in 
my view, the advent of realistic detailed simulations, in particular radiation (magneto-) hydrodynamic 
simulations casting Sc andinavians and Robert F. Stein as coauthors. Of cour se, I allude to the 



Carlsson & Steinl ( 1997) aco ustic p ropagation simulation, the ISkartlien et al.l (2000) acoustic excitation 



simulation, and the Ros enthal et al.1 (2002) shaken-fiuxtube and canopy-interaction MHD simulation. An 
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important strength of these simulations is that they are unusually close to observations, firstly in aiming to 
reproduce particular observed phenomena rather than prove grand presupposed concepts, secondly in the 
use of elaborate observation-like diagnostics to analyse simulation physics rather than just showing overall 
results. 

The issues discussed here suggest the following simulation to-do items: 

- shaken-fluxtube Ca II H & K emission; 

- cool-star Wilson-Bappu relation; 

- umbral flashes and running penumbral waves; 

- granular overshoot and gravity wave excitation; 

- gravity wave interaction with fluxtubes and canopies; 

- cool-star basal flux evaluation. 
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